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We discuss the effective donor/ acceptor coupling for a 
bridged electron transfer system [h| with a site-diagonal dis- 
order of bridge energies. The average spectral properties of 
the system are discussed by using the Wegner model (An- 
derson's type tight-binding Hamiltonian (TBH)) for the elec- 
tronic part of the problem. Spectral properties of the system 
are discussed using the concept of the functional inverse of 
the resolvent ("Blue's function", Q) for various limits of noise 
versus site-site coupling ratio. 



I. INTRODUCTION 

Electron transfer (ET) processes play a fundamental role 
in chemistry j^] and biology Q] with the range of their 
abundance ranging from photosynthesis and oxidative 
phosphorylation to molecular electronic design. In nu- 
merous biological examples of ET reaction, a single elec- 
tron is tunneling in an inhomogeneous medium over large 
distances of several angstroms. The intervening medium 
can be either a protein backbone or a sequence of co- 
factors embedded in a protein matrix. Due to a large 
separation between the donor and acceptor, direct elec- 
tronic coupling between the chromophores is negligible, 
rendering thus the question on the effect of medium on 
enhancement of the electronic coupling 0. The typi- 
cal suggestion is that the ET is mediated through the 
medium which is acting as a bridge which provides virtual 
states for the tunneling electron [||. In such a superex- 
change mechanism, effective electronic coupling between 
the initial and final states depends on the structure and 
flexibility of the bridging system. 

In this paper, we adopt a standard "effective two-level 
system" approach towards ET problem [§,0J|] . The tech- 
nique has been elaborated in the past years by formal- 
izing ET in the energy domain by use of diagonalization 
procedures which lead to transition amplitudes expressed 
in terms of Green function ]s|^,?]. To include possible 
fluctuations in bridge energies, we assume that the bridge 
Hamiltonian can be " sampled" from the class of random 
Hamiltonians. The underlying argument is that inhomo- 
geneous, protein medium can affect the site energies of 
the bridge and influence the kinetic ET rate. Our studies 
are performed in the limit of N — > oo where N stands for 



number of bridge orbitals. For such an extended system, 
we study effect of the noise on density of bridge electronic 
states. 

The analysis is performed within the framework of the 
random-matrix theory [ 1^] which turned out to be quite 
general and a powerful phcnomcnological approach to 
a description of various phe nomena such as quantum 
chaos [jll| , complex nuclei , chaotic scattering |l3| and 
mesoscopic physics jl4|] . Aspects of vastly different phys- 
ical situations such as electron localization phenomena 
in disordered conductors and semiconductors fUSl , disor- 
dered quantum wires |ll| and quantum Hall effect JTt| ] 
can be described in the language of the random matrix 
theory. In all the realms mentioned above, the Hamil- 
tonian of the system is rather intricate to be handled or 
simply unknown. In such cases the integration of the 
exact equations is replaced by the study of the joint dis- 
tribution function of the matrix elements of the Hamil- 
tonian P(H). 

The natural way of addressing the problems of random- 
ness coupled to various sources is to use technique of 
the "free random variables" IH-H. This method pro- 
vides an elegant way of "linearizing" the process of deter- 
mining the average eigenvalue distributions for convolu- 
tions representing an analogue of the logarithm of Fourier 
transformation of the usual convolutions. Recently, the 
generalization of the "addition law" for hermitean ran- 
dom matrices to the non-hermitean case has been de- 
rived |2l],|2^] and applied by us to study properties 
of a dissipative two-level system. Similar diagrammatic 
approach has been used by other authors to inves- 
tigate spectral properties of the Fokker-Planck operator 
that describes particles diffusing in a quenched random 
velocity field. 

The paper is organized as follows. In Sections 2 and 3 we 
discuss briefly the model Hamiltonian reporting results 
known from the literature. Spectral properties of the 
bridge Hamiltonian are further estimated in the large N 
limit. In Section 4 a model of disorder superimposed on a 
tight binding Hamiltonian of Sec. 2 is presented. By use of 
the concept of " free variables" , we estimate Green func- 
tion of the disordered system. Average properties of the 
system can be further inferred by studying the structure 
of the distribution of eigenvalues of the bridge Hamilto- 
nian which follows derivation of the Green function for 







the model. 
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II. PARTITIONING TECHNIQUE, EFFECTIVE 
COUPLING AND GREEN FUNCTION OF THE 
BRIDGE. 

The transferring system is assumed to have an electronic 
mrt described by a tight-binding (Hiickel) Hamiltonian 



(2.1) 



H = e(\D)(D\ + \A)(A\) + H coupl + H bridge 



where e stands for the energy of both donor and acceptor 
states l-D), \A) (for simplicity they are taken here to be 
equal) and the coupling H coup i and bridge Hb r id ge Hamil- 
tonians are given by 



iV 



H coupl = (pDiWibil+PAiWihl) + he (2.2) 



and 



N 

Hbridge = + / ] Pij \ h) (bj \ 

i ijtj 



(2.3) 



with 0ij being a symmetric site-diagonal matrix. 
The effective Hamiltonian is obtained by the standard 
Wigner-Weisskopf reduction |t],^6|,^7J by partition. In 
this method the total Hilbcrt space of the problem is di- 
vided into two subspaces (spanned by donor and acceptor 
states and bridge states, respectively) and the Hamil- 
tonian is integrated over the elements of one of them 
(here, the bridge states), eventually mapping the eigen- 
value problem of a high dimension onto a lower one. The 
time-independent Schrodingcr equation for the partition 
scheme is |E7f 



(H - ET)c = Mc = 
with matrix M given by 



M , M aa M ab 



M ba M bb 



(2.4) 



(2.5) 



By eliminating subspace "b" (spanned by intervening 
bridge states |i>i)), the effective Schrodinger equation be- 
comes 

(M aa - MoiM^MhOc = M QQ c a = (2.6) 

so that the reduced effective Hamiltonian of the system 
is the 2x2 matrix 



Heff = M Q 



Ela 



(2.7) 



with energy E set up (in zeroth order) to e, i.e. to 
donor/acceptor energy. The formalism yields the effec- 
tive coupling given by [pLpq] 



where bracketed expression stands for the Green function 
of the bridge 



G( 6 ) — (H br idge — e ) 

with H br idge given by the TBM, 

^bridge = e bSmn + P(dm,n+1 + S m ,n-l) ■ 



(2.9) 



(2.10) 



The Hamiltonian above describes a system with N local- 
ized sites, each with energy e b and a hopping between the 
neighboring sites with hopping parameter /?. 

Fr om t he spectral representation of the Green's func- 
tion (|2^) 



with eigenvalues 



e b + 2(3 cos ( 



N+T 



and spectral coefficients 

, 2 



) 1 / 2 sin( 



AT -Fl- 



it follows, that 



de 



K«0 

A - e 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



where in the limit of N — > oo we have replaced the sum 
in (2.9) by an integral with the level density v(e) as a 
weight. The localization length 7 -1 measuring the decay 
of an eigenvector ove r the chain composed of N units can 



be deduced from eq.(|2.11 



log/3 



and in the limit of an infinite chain reads 
l{tu)= I de v(e)\og\e v - e\ - log/3 



(2.15) 



(2.16) 



bringing the dependence on the density v(e). 
By changing to complex variable, the density of eigenval- 
ues can be conveniently defined in terms of the trace of 
the resolvent of the Hamiltonian 



G ( e ) = ^ Tr ( ^ 

1\ \ e — tl b ridge 



(2.17) 



The density of states for H br id ge is then given by 
1 

1/(6 



L Tr (S(e - H brldge )) = ~ImG(e + iX) (2.18) 



N IT 

and follows from the discontinuities of the resolvent 
( [2.17 ) along the e-axis. 
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III. KINETIC RATE. 

The kinetic rate for the electron transfer mediated by the 
bridge can be evaluated according to the definition |j],|]|5) 



Ht) = j t \{<l>f\Mt))\ 2 

where \faf) is a final state and 



(3.19) 



(3.20) 



with ipi sta nding for the exact eigenstatc of the Hamil- 
tonian (2^2). The latter can be decomposed into H = 
Hq + Hi in which Ho describes diagonal part of the total 
Hamiltonian and Hi defines interaction inducing charge 
transfer between the donor and acceptor ||. The exact 
eigenstates of H can be represented as 



ipi = 



1 



H + iX 
1 



Hii>i 



-Hie 



d-H + iX 

where fa is the initial state of the system and 

1 



T = Hi+Hi 



H + iX 



Hi 



(3.21) 



(3.22) 



stands for the transition operator which to the lowest 
order of approximation is given by 



T = Hi+H 



1 



1 €i - H + iX 



Hi 



(3.23) 



By use of the transition operator, the kinetic rate ( 3.19| ) 
can be rephrased in the form 



2tt 



(faflTlfajfSief-a) 



2tt 



(faf\Hi\^) r«(e/-ei). 



(3.24) 



By assuming that the vibrational relaxation in the molec- 
ular system is much faster than the electron transfer, the 
above equation can be expressed as 



k = irEE P « !<*/»' \T\\^v)\ 2 S(e fv 



,) (3.25) 



where (if) and (v, v') stand for electronic and vibronic 
states, respectively and Pi v is Boltzmann weight factor. 
In the adiabatic approximation, the wavefunctions of the 
system can be written as products of electronic and vi- 
brational wavefunctions 



(3.26) 



Substitution of the above approximation to the evalua- 
tion of matrix elements of T yields 



Q> fv »\T\iM = (Q fv „\H lfi \Q iv ) (3.27) 

(®f v "\Hlf m \Q mv r) (Q m v'\V m i\Qiv) 



E 



+ iX 



In the above formula, elements Hi f m represent elec- 
tronic matrix elements between states $/ and (fr m . Con- 
don approximation |^-|5j allows for factorization of these 
terms from the Franck- Condon factors (@fv"\®iv)i an d 
in the case when the electronic gap between the donor 
and acceptor states is larger than vibrational energies in 
the system, leads to the ET rate 



2tt, 

T 



k ^IV 2 ^ 



v' ,v 



Piv\(Qfv'\Qiv)\ S(e fv 



.,)• (3-28) 



The rate constant can be thus written as a product of 
electronic part and nuclear part, averaged over vibra- 
tional states. Throughout the paper we will focus on 
electronic deg rees o f freedom and estimation of the tran- 
sition matrix ( 3.27 ). 

Note, that by neglecting a direct electronic coupling be- 
tween the donor and acceptor and interpreting 



($ f \Hi\<l> N ) = (fa\Hi\$i) =/3i 



(3.29) 



the formula for the transition mat rix (3.27) is identical 
with the the effective coupling (2^8) and for the chain of 
N elements of the bridge reduces to ID 



T fi = H DA = (-1) 



a 2 R N-1 
N Pi P 



Jf - -filGiN 



(e 6 -e) 

where the element Gin of the Green's function is 
(-l) N + 1 p N ~ 1 



G 



IN 



detHf, 



■ridge 



(3.30) 



(3.31) 



with det sta nding for the determinant of the bridge 
Hamiltonian (2.10): 



detN — 
with 



2V 



a* 



(3.32) 



2ft 



(3.33) 



In the N — > oo limit det approaches ft N for \a\ < 1 and 
oo otherwise, thus the transfer element Gin is 1//3 for 
a small donor-bridge energy splitting (|a| < 1), and 
for tunneling energies outside the band of the width 2ft 
(the latter corresponds to the Bloch-like extended states 
of the bridge @). 

Localization constant introduced in the former para- 
graph is related to the decay of the transfer matrix |28[ 



7 = - Jim T7 m |Gi,jv| 

JV^oo JSI 



(3.34) 
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and, in the case of disordered bridge, would require usin g 
the density function v(e) to estimate the average ( 2.16 ). 
In the deterministic case (no site diagonal 'disorder') one 
has 







max log \a ± v? 



II lal > 1 



(3.35) 



so that the states with energies lying within the band 
(—2/3, 2(3) are localized over the infinite range. 



IV. TRANSITION THROUGH THE 
DISORDERED CHAIN. RANDOM MATRIX 
MODELS. 

As it stands, the bri dge Hamiltonian (tight-binding 
Hamiltonian TBH, eq. ( |2.3| )), can be generalized, as in 
the case of the Anderson model, to include disorder of the 
bridge chain. In what follows, we will adopt the model of 
the site diagonal disorder which is equivalent to the ones 
studied by Wegner (§§ and JVeu and Speicher 
The idea of Wegner was to generalize the Anderson model 
by putting n electronic states at each site of the d— 
dimensional lattice and describing the disorder by Gaus- 
sian random matrices in the electronic states. For n = 1 
Wegner's model reduces to the usual unsolvable Ander- 
son model and becomes exactly solvable for n — > ocjj]. 
Further generalization of Wegner's model has been dis- 
cussed by JVeu and Speicher who have used a mathemat- 
ical concept of "freeness" to use a more general ensem- 
bles of random matrices allowing arbitrarily distributed 
disorder. "Freeness" , introduced in mathematical liter- 

19|j3|] has 



ature by VbicuJescu, Pastur and Speicher 
been also popularized recently in physical applications by 
Brezin, Zee, Janik et al. 

The formalism can be translated to describe spectral 
properties of a Hamiltonian of the form 



H = H 



D 



H 1 



(4.36) 



where H D is deterministic, and H R random part of the 
operator. By assuming that H D and H R are free with 
respect to average over the disorder, VbicuJescu, Pastur, 
JVeu and Speicher have shown that the diagonal part of 
the one-particle Green function associated with the total 
Hamiltonian H satisfies equation 



G(e) =G D {e-Y,(G(e))} 



(4.37) 



where the argument of G D is e — £ with E being nothing 
but the self energy determined by 



1 The method described in that paper also applies to 
the original n — 1 Anderson model and then it is equivalent 
[po|j3l| to the CPA ("coherent potential") approximation 
which works remarkably well. That allows us to extrapolate 
the results to finite n configurations. 



G = 



£(G«) 



where 



G 1 



N 



-Tr 



-H R 



(4.38) 



(4.39) 



The same result has been rederived by Zee [£| who 
through his diagrammatic analysis, introduced the 
"Blue's function" that is just the functional inverse of 
the resolvent 



B[G(e)} = e 
and satisfies the additivity law 



B 



D+R 



(e)=B JJ (e)+B R (e)--. 

e 



(4.40) 



(4.41) 



Both equations ( 4.37 ) and ( 4.41 ) coincide if one identifies 
5(e) = E( £ ) + e -^ _ 

Our further analysis is based on the assumption that 
nodal energies of the bridge are randomly distributed 
with off-diagonal elements of matrix Hb r id ge being con- 
stant. In the limit of large JV, deterministic resolvent of 
the bridge Hamiltonian yields 



G 



D i 



1 N 

-y 

N ^ 



1 



N ' — ' £ 

fc=l 



2/3 cos ^ 



tlx 



2(3 COS 7TX y / e 2 _ 4/? 2 



(4.42) 



(note, that the average et — has been set up to 0) with 
the deterministic "Blue's function" 



B" = j\+4(3 



(4.43) 



Evaluation of the deterministic Gin element for this case 
leads, in the limit of (et — e)/2(3 > 1, to the usual Mc- 
Connel result 



lim Gin = . 



N- 



(4.44) 



Qualitative estimate of the transfer matrix for the "ran- 
dom plus de term inistic" case follows now, according to 
the formula (3.31), determination of the inverse determi- 



nant for the full (diagonal random plus TBH determin- 
istic) bridge Hamiltonian: 



/ {-1) N+1 P N - 1 
\det{H D + H R ) 



Gi 



JV 



(4.45) 



For (3 "large" det is dominated by (3, giving det = (3 N ; 
for small (3 the determinant is dominated by the random- 
ness of the diagonal elements in the bridge Hamiltonian. 
Thus naively, in the limiting case of "large" (3 (and for 
an infinite bridge) (3 in the numerator is expected to be 



3 



negligible compared to the det leading to Gin propor- 
tional to 1/(3. More careful analysis (see below) shows, 
that randomness always increases 7, leading to the faster 
decay of the transfer matrix Gin over the distance. 

If randomness of e& is chosen to be sampled with the 
semi-circle lawr 



£»(e&) 



I 



27TCT 2 



4<7 2 



(4.46) 



which is equivalent to a random matrix model |M , ^ , |9[ 
of the noisy contribution H R to the bridge Hamiltonian, 
Gre en's f unction associated with H R takes the form (see 
eq. (p7|)) 



G R {e) 



e =F Ve 2 - 4cr 2 



(4.47) 



.-<j 2 G R (e) 2a 2 
The functional inverse of G R is given by Blue's function 



B 



e 



(4.48) 



and, after using the addition law ( 4.41 ) leads to the fol- 
lowing equation for the Green's function of the system 



cr 4 G 4 - 2<xr 2 G i + (e 2 - Af3 2 )G 2 -1 = 0. 



(4.49) 



The end-points of the spectra may be calculated from the 
equation 



dG 
de ' 



(4.50) 



or, equivalently from the discriminant of the above equa- 
tion. By introducing the rescaled variables 



a 



/3G-> G 



(4.51) 



the discri mina nt (end-point condition) for the Pastur 
equation ( [l.49| ) is 



4e«+(<7 4 - 48)^ 



- 16(12- 5<7 2 )e 2 -16(4+cr 4 ) 2 : 



0. (4.52) 



In the limit a — ► that reduces to (e 2 — 4) 3 = 0, thus the 
support of the spectrum is one interval, with end-points 
±2, while in the limit a — > 00 ((3 — ► 0) one gets once again 
one interval with end-points e = ±2<r. Generally, the 
discriminant equation has only one pair of real solution 
for any value of cr. 

The spectra for different values of a is shown in Fig. [j]. 
'Deterministic' spectrum is divergent close to the the end- 
points of the band. For an increasing value of the noise 
intensity cr, the distribution v(e) flattens and, eventu- 
ally changes from the bimodal function with two peaks 
located at the ends of the support to the unimodal dis- 
tribution with a broad hump around the center of the 



We note that choosing other randomness does not change 
our results qualitatively. 



support. The curvature around zero changes sign for a 
critica l valu e of the noise intensity which can be found 
from (h4£) analyzed around e = 0. With the scal- 
ing ( 4.51 ) one gets then the critical value of er c in units 
of & 



12 



1/4 



1.8612. 



(4.53) 
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FIG. 1. Eigenvalue distribution for the TBH of the bridge 
for different values of the noise/deterministic ratio, cr. 

Note that from the characteristics of v[e) it follows, 
that the localization length j^ 1 has a discontinuity at 
the end-points of the spectra for a = 0. For energies 
within the band, 7 is zero and switches to a finite value 
after the condition \a\ > 1 is reached. The noise in some 
sense softens the abrupt change of the 7 at the endpoints 
of the support (c/.Fig.2) and shortens the localization 
length within the band. 

The "deterministic plus random" Gin can be evalu- 
ated (in the limit of large N) following Neu and Spe- 
icher fill, 



Gw(e) = Gf JV (e-S[G(e)]) 



(4.54) 



i.e., the same Pastur equation ( |4.37| ) holds also for the 
non-diagonal elements of the t otal Green's function. The 
general analysis of Eq ( 4.54 ) shows, that the solution 
7 = can not be achieved any more for any combination 
of £{, and (3. However for "small" disorder the states with 
energies (ej — e)/2/3 < 1 are stable with small decay coef- 
ficient (see Fig. g). The most stable energies are sampled 
around e = with decay coefficients 



7 : 



log -L f yj ^l-h7 4 /2 - 1 + ^1^/2 + 1) ( 4 . V, ) 



yielding 7 = only for vanishing randomness (cr = 0). 

In Fig. U one may note the cusp-like discontinuities 
around the end-points e e of the spectrum. Their appear- 
ance may be understood analytically in the large cr limit, 
where by evaluating the resolvent G we neglect the de- 
terministic part of the Hamiltonian. The resolvent then 
reads 



G%) = 



4ct 2 



2a 2 



(4.56) 
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with end-p oints e e = ±2(7. Evaluating (4.54) and insert- 
ing it into ( 3.34 ) the expansion around the end-points e e 
for e > e e leads to 



7 « log 



2 

showing the observed "cusp" around the end-point 

1.5 



(4.57) 




FIG. 2. Decay coefficient 7 as a function of energy for var- 
ious intensities of the noise. (3 has been set up to 1. The 
completely deterministic system (a — 0) has 7 = for |e| < 2. 



V. CONCLUSIONS 



matrix Tu. Thus the diagonal noise reinforces the ex- 
ponential decay of the effective coupling with increasing 
distance between the donor and acceptor. 
The method of Blue's function applied here is also 
suitable for nonhcrmitean ensembles of random ma- 
trices which are used in quantum theory of dissipa- 
tion (Tins Thc 

paper discusses so far only the case 
of the hermitean ensemble. Extension of the formalism 
to models of dissipative transport will be presented else- 
where. 
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In this paper we have discussed briefly properties of a 
bridged electron transfer affected by site diagonal dis- 
order in the bridge Hamiltonian. Instead of using the 
standard language of the Green's function formalism, we 
have chosen to work with the functional inverse of the 
Green's function ("Blue's function"). The method stems 
from the powerful mathematical concept of free random 
variables (ijjnj which has been shown to be an elegant 
tool in various applications of the random matrix the- 

ory|HI§. 

Our model refers to the situation when the bridge is 

long enough to be considered infinite. Site diagonal disor- 
der is assumed in the form of a random matrix model re- 
sulting in placing at each site of the bridging chain a ran- 
dom matrix with a semi-circular distribution of energies 
(Wegner model, p{|). That would correspond to a situ- 
ation where the bridge energies are n-electron functions. 
Presence of noise extends and flattens the spectrum of the 
TBH Hamiltonian. For donor/acceptor energies sampled 
from the center of the TBH band, the noise increases the 
inverse localization length 7 and leads to a fast decay 
of the electronic coupling with the distance measured in 
bridge units. The kinetic rate becomes thus vastly re- 
duced by the noise (note that our discussion relates only 
to the electronic part of the ET rate). The distribution of 
TBH energies has been shown to posses a noise-induced 
characteristics which depends on the value of the critical 
noise intensity a c . As expected p9| , p2 34 1, the diagonal 
disorder localizes the eigenfunctions of thc TBH Hamilto- 
nian resulting in the reduction of the electronic transfer 
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